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CN . Abstract 
^ ■ 

Ph! a La grangian is introduced which includes the couphng between magnetic moments m and the 



degrees of freedom <t of a reservoir. In case the system-reservoir couphng breaks the time reversal 
symmetry the magnetic moments perform a damped precession around an effective field which is 



. self-organized by the mutual interaction of the moments. The resulting evolution equation has the 



form of the Landau-Lifshitz-Gilbert equation. In case the bath variables are constant vector fields 
the moments m fulfill the reversible Landau-Lifshitz equation. Applying Noether's theorem we find 
■ conserved quantities under rotation in space and within the configuration space of the moments. 
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I. INTRODUCTION 



The dynamics of magnetic systems is described in a wide range of time and length scales 
from a quantum approach up to a macroscopic thermodynamic access. On a coarse-grained 
mesoscopic level the relevant electronic degrees of freedom are grouped into effective mag- 
netic moments. As the consequence the magnetization is characterized by a spatiotemporal 
vector field m(r,t). Owing to the mutual interaction between the magnetic moments they 
perform a precession motion around a local effective field giving rise to the propagation of 
spin-wave excitations. Due to a system-reservoir coupling the precession of the moments 
should be a damped one. To analyze this situation one has to specify the coupling between 
the system and the bath. The most popular approach to incorporate dissipation is the em- 
bedding of the relevant system into a quantum-statistical environment which is assumed to 
remain in thermal equilibrium. The reservoir is often represented by harmonic oscillators or 
spin moments which are analyzed by path integral techniques [l- 



. A specification of the 

path integral approach for spin systems can be found in [443]. A more generic description 
concerning dissipative semiclassical dynamics is presented in js]. Although the application 
of path integrals can be considered as an intuitive formalism analytical calculations are often 
impossible and numerical schemes are necessary. In the present paper we propose an alter- 
native way to include dissipative effects for mesoscopic magnetic systems. On this level the 
analysis of magnetodynamics is performed properly by applying the Landau-Lifshitz-Gilbert 
equation designated as LLG j^,llo|. A comprehensive survey of magnetization dynamics is 
given in . Our investigation can be grouped in the effort to understand dissipative mech- 
anisms in magnets. So, a non-linear dissipative model for magnetic systems was discussed in 
. On the relation between fluctuation-dissipation theorems and damping terms like that 



one occurring in the LLG was reported in 



131]. The dynamical response of ferromagnetic 



shape memory alloy actuators can be modeled by means of a dissipative Euler-Lagrange 



equation as performed in 



14l |. Likewise, the pinning of magnetic domain walls in multifer- 



roics is discussed in terms of the EL equations in 15|. An alternative ansatz is introduced in 

n 

[Iq . where a Lagrangian density is obtained based on a projection onto the complex plane. 
The procedure gives rise to a dynamical equation which is equivalent to the Landau-Lifshitz 
equation. Different to the mentioned approaches the present paper is aimed to derive an 
equation of motion for a magnetic system on a mesoscopic scale under the influence of a 



bath which hkewise consists of mesoscopic moments. Following this idea we propose a La- 
grangian comprising both fields, m(r, t) as the system variables and cr as the bath variables. 
The bath becomes dynamically active by the coupling to the system. In case the coupling 
between system and reservoir breaks the time reversal symmetry the motion of the moments 
m(r, t) is damped. The Lagrangian is modified in such a manner that dissipation can occur. 

II. THE LAGRANGIAN 

As indicated we are interested to construct a Lagrangian describing the motion of a magne- 
tization vector field within a bath of spins. This reservoir should infiuence the measurable 
magnetization due to the mutual interaction. Let us formulate the general assumptions for 
the underlying model. The magnetic order is originated by single magnetic atoms which 
occupy equivalent crystal positions. Here we refer to a continuous description in terms of 
a field vector denoted as m(r,t). Because the ferromagnet is considered below the Curie 
temperature a sufficient number of microscopic spins preferring a parallel alignment are in- 
cluded in m, i.e. the effective magnetic moment is given by m(r, t) = ^ ■ /Xj where the 
sum is extended over all microscopic moments within a small volume around the spatial 
coordinate r at time t. As each axial vector the moment fulfills m(— t) = — m(t). The bath 
in which the moments are embedded consists likewise of mesoscopic spins. They are denoted 
as cr and are also composed of microscopic moments 77 j. This bath moments which play 
the role of 'virtual' moments are also axial vectors changing their sign by time inversion. 
A further new aspect is that the coupling between the real and the virtual moments is not 
assumed to be weak. As the result the complete system consists of two subsystems. One of 
them abbreviated as LI is occupied exclusively by the real spins with the moments m and 
the other one denoted as L2 is occupied by the bath spins cr. The situation is illustrated in 
FIG. [H Now let us introduce the action 



where the set {qa} consists of the set of both moments cr and m. The Lagrange density 
comprises three terms 




(1) 



£[m, m, Vm 



• (T, Vo-] = [m, m, Vm] + [Vct] + C^"""^ [cr, m] , 



(2) 



3 



jfXT] jr]r] jufi 




Vj-k 'Hj fJ'i Vj+l Vj+l+k l^i+l 

FIG. 1. (Color online) Schematic illustration of the basic model. The red spins represent the 
magnetic moments /ij and refer to the lattice LI introduced in the text. The green spin vectors 
r]i build the bath lattice L2. Interactions are possible between the /^j and rj^, /ij and /Xj and ?7j 
and r]y The respective coupling strengths correspond to the coupling parameters in Eqs. ([3])- ([5]) 
as follows: J^''^ ^ J^'^' o and J™ o J^-^). 



where Z^^™-) indicates the Lagrangian of the magnetic system, C'^"^ represents the reservoir 
and the interaction term is denoted as To be more specific the magnetic moments of 

the system interact via exchange couphng defined by the Lagrangian 

where J^^^^ is the couphng parameter, diagonal in the isotropic case. The first term represents 
the energy density of the magnetic system. Because we are not considering the acceleration 
of magnetic moments a term of the order rh^ is missing. Moreover, the magnetic moments 
perform a precession around an effective magnetic field, which is self-organized by the mutual 
interaction. Therefore the vector potential A depends on the moments, i. e. A = A(m(r, t)). 
The coupling has the same form as the minimal coupling in electrodynamics. The bath 
Lagrangian is defined in a similar manner as 

with the coupling constant J^^p ■ Eventually, the interacting part between system and bath 
is written as 



4 



with the couphng strength ■ The second term is constructed in the same manner as 
in Eq. (|3]), where the potential B(cr) will be specified below, see Eq. ([8]). The dynamics 
of the bath variable <t remains unspecified for the present, i.e. the Lagrangian does not 
include a term of the form oc <t. Owing to the constraint, introduced in the next section, 
the dynamically passive bath is sensitive to a change of the system variables m in such a 
manner that small variations of the system variables m are related to small variations of 
cr. This procedure leads to a coupling between bath and system so that the time reversal 
symmetry is broken. 



III. RELATION TO THE LANDAU-LIFSHITZ-GILBERT EQUATION 

In this section we find the equation of motion for the magnetization m(r, t) from Eq. ([2]) 
combined with Eqs. (l3])-(l5l). Using the principle of least action it follows 



dL d dC 



da 







\axa J 



Sas + 



dC d dC d dC 



drrijs dt dirnp dx^ g ( 

\ dXa 



drrip = , 



(6) 



where and 6a ^ are the small variations which drive the value for the action out of the 
stationary state. In general, one derives a system of coupled partial differential equations. 
However, to proceed further let us impose a constraint on the system. A small variation of 
(T/3 should be related to a small variation of m/3. Thus, we make the ansatz 

6ai3 = —K 617113 , with K = const > . (7) 

Notice that this condition should be valid only locally but not globally. Insofar Eq. ([7]) is 
comparable to an anholonom condition in mechanics. Moreover relation ([7]) is in accordance 
with the behavior of the moments m and cr under time inversion. Physically the last relation 
means that the bath reacts to a change of the system only temporarily. Because the system- 
reservoir coupling should typically break the time reversal symmetry the expansion of the 
function By{a) in terms of a includes only odd terms. In lowest order we get from Eq. (jS]) 

Bi^{a) = —ca^, with c = const . (8) 

Due to Eqs. ([7]) and ([8]) the second term in Eq. ([5]) is of the form oc cr • m. Such a term is not 
invariant under time reversal symmetry t — )■ — t. As demonstrated below the broken time 



inversion invariance gives rise to damping effects. Inserting Eqs. ([7]) and ([8]) into Eq. ([2]) 
and performing the variation according to Eq. ([6]) we get 

= — + cap + Kcm^ 

- [ J^'") - K J("") ] V^m^ - [ J^™") - K J(") ] V . 

Here we have assumed for simphcity that all coupling tensors J are diagonal: Jajs = JSa/s- 
The first term on the right hand side in Eq. ([9]) reminds of the field strength tensor in 
electrodynamics [l^]. Thus, we rewrite 

OA, dAp' 



rhy = Fr^ Thy 



mx(VmXA(m)) . (10) 



As mentioned above the vector function A(m) is regarded as vector potential which depends 
on space-time coordinates via the magnetic moment m(r,t). In vector notation the last 
equation reads 

m X ( V™ X A) = [ J^'") - ft J("^") ] V'm - /tc m - c o" + [ J™" - « J^'^ ] V V . (11) 

If one is interested in weak excited states of a ferromagnet it is reasonable to assume that 
the direction of the magnetization in space changes slowly while its absolute value is fixed, 
that is = 1. Without loss of generality we have set the amplitude of m to unity. In order 
to proceed it is necessary to specify the condition which should be fulfilled by the function 
A(m). Having in mind the LLG then we make the ansatz Vm x A(m) = gm, g = const. 
Based on these assumptions we get from Eq. ( ITT!) 

dm 1 / -^-^ \ Kc I dm , . 
^ = -(^mxH„.j--(^mx_|. (12) 

Here the effective field is given by the expression 

da 



H 



eff 



Eq. f[T2|) is nothing else than the Gilbert equation 



10| by relating the prefactors as follows 



1 '^'^ f-iA\ 

7 = — , a = = Kcy , (14) 

9 9 

where 7 and a are the gyromagnetic ratio and the Gilbert damping parameter, respectively. 
Since both parameters are positive quantities it follows that g < as well as kc > 0. Further, 



Eq. (fT2|) can be converted into the form of the equivalent and widely used Landau-Lifshitz- 
Gilbert equation which reads 



dm 7 / TT N «7 

fm X Hcff) 



m X fm X H, 



cS 



(15) 



dt (l + «)2' (l + «2) 

Both quantities 7 and a are still related to the model parameters by the expressions in 
Eq. f ll4p whereas the the effective field Hes is given by Eq. f ll3p . Now we want to analyze 
this expression and in particular, to assign a physical meaning to the more or less ad hoc 
introduced quantity cr. In doing so one can distinguish four different cases: 
(i) The bath is not included which corresponds formally to cr is a constant vector depend- 
ing neither on coordinates nor on time. Then obviously all derivatives with respect to the 
coordinates and the time of cr disappear in Eq. (IT^ and consequently, the set of {qa} in 
Eq. dl]) does not include cr. From here we conclude that the variation fulfills 5cr = in 
Eq. which can be easily reahzed setting k = 0, cf. Eq. ([7]). Thus, the effective field 
in Eq. (fT3|) comprises the pure exchange interaction J^™'^ between the magnetic moments 
and the damping term in Eq. f|T2l) is absent due to a = in Eq. f|T4l) . A constant bath 
field cr lead to the Landau-Lifshitz equation in the exchange interaction approach without 



damping, compare 18|. It describes the precession of magnetic moments of an effective field 
which is self-organized by the mutual interaction of the moments. 

(ii) cr = cr (if:) depends only on the time and not on the spatial coordinates. Regarding 
Eq. f|T3l) the effective field is modified by two additional contributions, namely one propor- 
tional to V^m, originated in the exchange interaction of the magnetic moments, and the 
other one oc &. The latter one could be associated with an external time dependent field 
or, if cr points into a fixed direction, gives rise to magnetic anisotropy. In that case the 
anisotropy axis is spatially constant but the amount of the anisotropy is changing in time. 
Such a situation could be realized for instance when the ferromagnetic sample is excited by 
the irradiation with electromagnetic waves. As already mentioned the exchange coupling 
J^"^^ is supplemented by a term -fi:j("'^). In this manner the exchange interaction is influ- 
enced by the coupling between m and cr although the spatial dependence of cr is not taken 
into account explicitly. 

(iii) cr = cr{r) depends only on the spatial coordinates and not on the time. In this case we 
first recognize that the coupling strength J*^™-' in the term oc V^m is influenced in the same 
manner as in case when cr = cr(t), see the previous point. Different to the former cases the 



expression oc V^cr becomes important for the effective field in Eq. (fT3|) . The appearance of 
this term suggests that spatial inhomogeneities of the surroundings of the magnetic system 
represented by m have to be incorporated into the effective field. It seems to be reasonable 
that the origin of this term is an inherent one and should not be led back to external fields. 
As possible sources we have in mind local varying fields like inner and outer demagnetization 
fields as well as accessible fields created for instance by different local temperatures, 
(iv) (7 = cr(r, t) is the most general case. Then external as well as internal fields are captured 
in the model. Thus, the effective field in Eq. (fT3|) can be rewritten as 

Heff(r,t) = Hexch(r) + h(r,t), (16) 

where Heg consists of two parts. The term Hexch = {J^"^'^ — /tJ*'™'^^)V^m is due to the 
exchange interaction between the magnetic moments whereas h(r, t) represents other possi- 
ble influences as discussed under the points (ii) and (iii). The function h is related to the 
quantity cr by 

h(r, t) = -c + [ Z'"'^) - «: Z'^)] VV(r, t) . (17) 

Remark that the formerly introduced quantity cr is related to the physically relevant effec- 
tive field by the first derivation with respect to the time and the second derivation with 
respect to the spacial coordinates via Eq. (11 7p . This equation is an inhomogeneous diffusion 
equation which can be generally solved by means of the expansion into Fourier series and the 
assumption of accurate initial and boundary conditions which depend on the actual physical 
problem. 



IV. SYMMETRY AND CONSERVATION 



After regarding the special example of the LLG we proceed with the investigation of more 
general aspects. The Lagrangian density allows to discuss the behavior under space-time 
dependent group transformation. For this purpose we apply Noether's theorem [l9| to our 
model. To be more precise we consider the conservation equation 



20| 



d 



dC 



/3 



dC ■ 



0, 



(18) 



Here, the expression in the square brackets are the components of the Noether current /□,. 
The term d/dXa in front of should be interpreted as an implicit derivative with respect 



8 



to time and three spatial coordinates. The symmetry operations AA'q and A'^a will be 
specified below. With regard to the Lagrangian in Eq. ([2]) we introduce the components 
= {mx,my,mz,<7x,<Jy,<Jz) and their partial derivatives with respect to the independent 
variables di^^^a = d'^a/dXp. Since we examine an Euclidean field theory a distinction 
between upper and lower indices is not necessary. Eq. (ITS]) can be rewritten by using 
Eq. ([6]). This yields 

In this equation we distinguish between the time and space variables t and Xa explicitly. 
Eq. ( |T9|) is the basis for the application of the following symmetry operations. Now we study 
the rotation around a certain cLXlS clS cl relevant one. Here we select for instance the z-axis. 
Performing a rotation in coordinate space with the infinitesimal angle AB the change of the 
X and y-coordinates obeys 




At = , Axa, = AR^fsXp , AR= \ . (20) 



In the same manner one can perform the rotation in the configuration space of the moments 
m and cr symbolized by the before introduced vector = {m^, mj^, m^, cr^:, o"j^, a^}. The 
transformation reads A\1/q = AS'a/3(A$)\I'^, where the rotation matrix is a 6 x 6-matrix 
determined by the rotation angle A$. Because both rotations in coordinate space and 
configuration space, respectively, are in general independent from each other we find two 
conserved quantities. Using Eq. (|T9l) it results 



D^C = 0, T,C = 0. (21) 
Here the two operators and f ^ are expressed by 

d 



f _ c(m) _|_ C(o-) 



a 



(22) 



The quantity is the generator of an infinitesimal rotation around the z-axis in the coor- 
dinate space 

Lz = y-^ - x-^ , (23) 
ox oy 



9 



and therefore, it is identical with the angular momentum operator. The other quantities 
Si^^ and Si'^^ are the corresponding generators in the configuration space of the moments. 
They are defined as 

cM _^ 

onix only , , 

y ^24) 



^ ^da^ ""day' 

These operators reflect the invariance of the total magnetic moment m + cr under rotation. 
Moreover the system is invariant under the combined transformation expressed by and 
f^;, where Dz offers due to the coupling between system and bath variables as well as the 
breaking of time reversal invariance a coupling between magnetic moments and the angular 
momentum. 



V. CONCLUSION 



In this paper we have presented an approach for a mesoscopic magnetic system with dissi- 
pation. The Lagrangian consists of two interacting subsystems characterized by the active 
magnetic moments of the system m and the dynamically inactive moments of the bath 
denoted as cr. Both systems are in contact so that a small local alteration of the system 
variables m is related as well to a small change of the bath variables a and vice versa. 
Due to this constraint we are able to describe the system by a common Lagrangian which 
incorporates both degrees of freedom and their coupling. In case the bath variables are 
constant then the coupling between both systems is absent and the whole system decays 
into two independent subsystems. The magnetic moments m perform a precession around 
an effective field which is self-organized by the mutual interaction of the moments. If the 
coupling between both subsystems breaks the time reversal symmetry the related evolution 
equation of the moments m is associated with the Landau-Lifshitz-Gilbert equation which 
describes both the precession of magnetic moments as well as their damping. It turned 
out that the bath variable cr can be linked to the effective magnetic field which drives the 
motion of the magnetic moments. As consequence the motion of the moments is infiuenced 
by the additional bath degrees of freedom. This influence is formulated mathematically and 
is described by an inhomogeneous diffusion equation. Finally, we have found conservation 
laws by means of symmetry considerations based on Noether's theorem. Aside from the 

10 



expected symmetry transformation in the coordinate space and the configuration space of 
the moments, the analysis offers in a non-relativistic Euchdean field theory an unexpected 
coupling between both. This point deserves further consideration. Our approach could be 
also considered as starting point for a further analysis in magnetic and multiferroic systems. 
Especially, we are interested in more refined models which include for instance higher order 
couplings or anisotropy in the Lagrangian. In multiferroic systems one could study the case 
that the magnetic and the polar subsystem have their own reservoirs. 

One of us (T.B.) is grateful to the Research Network 'Nanostructured Materials' , which 
is supported by the Saxony- Anhalt State, Germany. 
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